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We present a theoretical analysis of the temperature-magnetic field-concentration phase diagram 
of the multiferroic Mni_^M^W04 (M=Fe, Zn, Mg), which exhibits three ordered phases, with 
coUinear and non-coUinear incommensurate and with a commensurate magnetic order. The middle 
phase is also ferroelectric. The analysis uses a semi-phenomenological Landau theory, based on a 
Heisenberg Hamiltonian with a single-ion anisotropy. With a small number of adjustable parameters, 
the Landau theory gives an excellent fit to all three transition lines, as well as the magnetic and 
the ferroelectric order parameters. The fit of the magnetic and ferroelectric order parameters is 
further improved by including the effect of fluctuations near the transitions. We demonstrate the 
highly frustrated nature of these materials and suggest a simple explanation for the dramatic effects 
of doping with difi'erent magnetic ions at the Mn sites. The model enables an examination of 
different sets of exchange couplings that were proposed by a number of groups. Small discrepancies 
are probably a consequence of small errors in the experimental magnetic parameters. In addition, 
using the Ginzburg criterion we estimate the temperature range in which fluctuations of the order 
parameters become important. 

PACS numbers: 75.25.-l-z, 75.10.Jm, 77.80.Bh, 75.80.-fq, 75.40.Cx 

Keywords: Multiferroics; Magnetism; Ferroelectricity; Magnetoelectric effect; Landau theory ; Critical phe- 
nomena. 



O 

o 



oo 
<^ 
o 



X: 



I. INTRODUCTION 

Type II Magnetoelectric multiferroics are materials 
which exhibit coexistence between certain types of long- 
range magnetic order and a ferroelectric order. These 
materials are usually characterized by a strong magneto- 
electric coupling between their electric and magnetic de- 
grees of freedom. The magnetoelectric effect enables the 
control of the electric polarization by a magnetic field, or 
the control of the magnetization by an electric field. The 
study of magnetoelectric multiferroics is thus of great 
interest in condensed matter physics, both from basic re- 
search and technological applications points of viewi^r— 
In recent years, the interest in this field has grown af- 
ter the discovery of new materials with a large magneto- 
electric effect, such as ThMnOg^ TbMnaOg^ NigVaOgr^ 
CuFe02,— and CoCr204.— In those oxides, ferroelectricity 
appears in conjunction with a noncollinear spiral mag- 
netic phase, which breaks spatial inversion symmetry, 
and therefore allows the appearance of an electric po- 
larization. 

There are two different approaches to the theoretical 
treatment of such noncollinear magnetoelectric multifer- 
roics. One approach is based on first principles cal- 
culations using the density functional theory (DFT)J^ 
The second approach constructs a model Hamiltonian 
dictated by symmetry considcrationsiiiii^ Different 
mechanisms for the magnetoelectric coupling can then 
be suggested ii^"— In this paper we develop a semi- 
phenomenological model for describing the magnetic 
phase transitions of Mni_^M^W04 (M=Fe, Zn, Mg) and 
the induced ferroelectric polarization. The model is semi- 
phenomenological in the sense that some of the param- 



eters can be deduced from existing experimental data, 
while the others arc purely phcnomcnological. The mul- 
tiferroic MnW04 is a natural choice for such an approach, 
due to the vast experimental data that exists in the lit- 
erature. 

MnW04 crystallizes in the wolframite structure, which 
belongs to the monoclinic space group P2/c with (3 « 
91°. The unit cell includes two magnetic Mn^+ ions with 
spin S ~ 5/2 and orbital angular momentum L = at 
positions ~ (0.5, y, 0.25) and T2 = (0.5, 1 — y, 0.75) (in 
units of the primitive lattice vectors) with y — 0.685J^ 
In zero magnetic field, MnW04 undergoes three succes- 
sive phase transitions at temperatures T^g sa 13.5K, 
Tjv2 ~ 12.3 - 12.7K and Tjvi ~ 7 - 8K to phases 
which are called AF3, AF2, and AFl, respectively^^"— 
According to neutron diffraction experiments,— AF3 is 
an incommensurate (IC) antifcrromagnctic phase with a 
coUincar sinusoidal structure, AF2 is an incommensurate 
antiferromagnetic phase with an elliptical-spiral struc- 
ture, and AFl is a commensurate (C) antiferromagnetic 
phase with a coUinear ftii structure. The propagation 
vectors are Qjq = (—0.214, 0.5, 0.457) (in units of the 
primitive reciprocal lattice vectors) for AF2 and AF3, 
and qpi 2 = (±0.25,0.5,0.5) for AFl. In AF3 and AFl, 
the magnetic moments of the Mn^"*" ions align along the 
easy axis of magnetization, which lies in the ac-plane and 
forms an angle of « 35° — 37° with the a axis. Different 
studiesiiii^ reveal that a ferroelectric polarization, which 
is oriented along the b axis, develops in the AF2 phase. 

As opposed to MnW04, other isomorphic wolframite 
structures like FeW04, C0WO4 and NiW04 show only a 
single magnetic phase transition to a simple commensu- 
rate antiferromagnetic phase with the propagation vector 
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q = (0.5,0, 0)ji^ Those observations suggest that unUkc 
the isomorphic structures, MnW04 constitutes a highly 
frustrated system with complex competing interactions. 
The competition between the different interactions man- 
ifests itself in the sensitivity of the phase diagram to dop- 
ing with different transition metal ions at the Mn sites. 
It turns out that a small Fe concentration suppresses 
the ferroelectric phase AF2 and expands the stabilization 
range of AF3 and AFli^"— In contrast to Fe doping, it 
has been reported^^ that a small Co concentration sta- 
bilizes the ferroelectric phase at the expense of the AFl 
phase. A quantitative and microscopic understanding of 
the effect of Fe and Co doping on the multiferroic prop- 
erties and the phase diagram of MnW04 is quite com- 
plicated, since the exchange couplings of the M-M and 
M-Mn (M=Fe, Co) interactions, as well as the anisotropy 
parameters are not known. In order to overcome some 
of these problems, a much simpler magnetic system has 
been achieved by the partial substitution of Mn ions by 
the non-magnetic ions Zn^"*" and Mg^+i^i^ Those stud- 
ies reveal that the AFl phase is strongly suppressed as 
a result of magnetic ions dilution by non-magnetic sub- 
stituents. 

The frustrated nature of MnW04 was demonstrated 
by Ehrenberg et al.^ Using inelastic neutron scatter- 
ing they extracted 9 exchange couplings Ji — Jg for the 
superexchange interactions among the Mn ions. Later, 
Tian et al^ proposed different values for the 9 exchange 
couplings based on DFT calculations. Those values de- 
pend on an unknown on-site repulsion energy. More- 
over, the authors have noted that generally DFT calcu- 
lations tend to overestimate the magnitude of exchange 
interactions 3i Recently, the experimental data have been 
expandedi^ In that study. Ye et al. suggested some cor- 
rections for the values of the exchange couplings, and in- 
cluded two additional ones, J^g and J^. The two sets of 



experimental exchange couplings are summarized in Ta- 
ble |T1 The model we describe may help to compare these 
different sets of exchange couplings, by examining their 
consistency with different experimental observations. 

The outline of the paper is as follows: in Sec. |lT] we 
define the model. In Sec. Illll the results of the model are 
derived. In Sec. IIVI the model parameters are fitted by 
comparing its results with different experimental obser- 
vations. Here we compare the two sets of experimental 
exchange couplings with the fitted parameters. In Sec. 
IVl the Ginzburg criterion is applied to the specific case 
of the multiferroic MnW04, in order to examine whether 
the mean-field theory approach is valid. We conclude in 
Sec. IVII with a brief summary. 



II. THE MODEL 

In this section we develop the semi-phenomenological 
model. The spin Hamiltonian consists of a Heisenberg 
term with a single-ion anisotropy, which favors an easy 
axis in the ac-plane. According to experiments, the spin 
component along the hard axis in the ac-plane does not 
order in any of the phases. Furthermore, the transi- 
tions are almost not influenced by an external magnetic 
field along the hard axis. Hence we omit the hard axis 
component from the calculations and write the spin as 
S{R + t) = S^{R + t)x + S^{R + r)b, where x de- 
notes the easy axis in the ac-plane and b denotes the 
axis perpendicular to the ac-plane. Here S{R + t) is the 
thermal average of the dimensionless classical spin at po- 
sition R + T, where i? is a lattice vector and r is one of 
the two basis vectors Ti, T2 in the unit cell, indicating 
the locations of the Mn^+ ions. We study the following 
Hamiltonian: 



J 
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RRT-r'^T^ 



J{R+t,R! + T')S{R- 



S{Ii+r')-\DY^ 

R T=Ti,r, 



S,\R + t). 



(1) 
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Here J{R + t,R! + t') is the superexchange interaction 
energy which couples the spins R + r and R! -\- t' , 
and £) is a positive single-ion anisotropy energy. To find 
an expression for the magnetic free energy of the system, 
we expand the entropy in the spin components up to the 
fourth order 

TS=-\aTY, E S\R + t) 
R -r^T^.T^ 

-&E E S\R + t), (2) 
R T=ri,T, 



where a and h are positive parameters, and T is the tem- 
perature. Equation ([2]) gives the entropy relative to the 
high temperature paramagnetic phase (denoted by P) 
and thus the expression is negative. Combining Eqs. ([T]) 
and ^ we obtain the magnetic free energy 
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TABLE I: Superexchange couplings for the Mn^^ ion at = (0.5, 0.25) according to dilTerent inelastic neutron scattering 
studies. We denote z — I — yjW — 2 — y and u — 1 + y. The values are presented in units of kgK^ 
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F = - 

mag rj 



E E E x:^(H + T,i?' + r')5„(i? + x)5^(H' + T') + 65^ 



J2 S'{R + t), (3) 



where the 4x4 inverse susceptibility matrix is block di- 
agonal 

X^liR + T,R' + t') = [{aT - DJSj^ j^,S,,,, 

-J{R + t,R' + t')]S,^P, (4) 

with Di = = D and D2 = D^^ = 0. Below, we exploit 
the Fourier transforms of the spin components, 

R 

S^iR + T)=Y^ S^{q, x)e-9-(^+^) . (5) 
Q 



Here q is in the first Brillouin zone and N is the num- 
ber of unit cells. In terms of the Fourier transform, the 
magnetic free energy per unit cell, f = F /N, is: 



/™a,= 2 E E E^^;^(9;r,r')5„*(q,r)5^(q,r') 

T, T'=Tj. X2 Q,;3 = l q 

+ ^E E E e-^G-^[S,iq„T)S,iq„T)S,{q„T)S,{q^,T) 

G T=Ti,T^q„q2,q^,q^ 

+ 5'fc(gi, -r)5fc(q2, -r)5'b(q3, ^)Sbi(l4, t) + 25^(qi, r)S'^(q2, T)S'f,(q'3, T)5'b(q4, r)] S^q^ + + Q3 + Qa-G), (6) 

I 

where G is a reciprocal lattice vector and the Fourier In the last expression the sum is over all lattice vectors 
transform of the inverse susceptibility matrix is given by R. The four eigenvalues of the matrix ([7]) are 
the block diagonal hermitian matrix 

X-0{q;T,T')=[(aT-DjSr^r' 

- J(q;r,r')]<5„,^, (7) 

with J(<7; T, t') being the Fourier transform of the 2x2 
matrix J{R + t,R! + t') 



J{q- T, r') = ^('^' R + T')e-"l<R+'^'~'^^ . (8) 
R 



C±Jq,T) = ciT-D-X^{q), 
C±,b(q,r)=ar-A±(g), (9) 
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and the corresponding eigenvectors are 



1 
71 



1 

71 



V 



(10) 



Here, are the two eigenvalues of the matrix (|8]) and 
4>{q) is the phase of J(q;l,2). Assuming 11 exchange 
couphngs as in Ref. [H, these two eigenvalues are given 

by 

A±(g) = ±2^K^\q) + K^^iq) + 2 cos{2^qb)K^{q)K.^{q) 
+ 2A,{q), (11) 
with the following definitions: 

Ai(q) = J3 cos(27rq^) + J4 cos(27r(3'^) 

+ J5 cos(27rgf,) + cos [27r(g^ + gj] 
+ cos [27r(q„ - gj] , 

A2(q) = Ji cos(^<zJ + Jg cos [2^(q, + ^)] 



+ J7COS [27r(g,-|)], 
^3(9) = J2 cosinqj + Jg cos [27r(q<, + |^)] 
+ J9COS [27r(g,-|)]. 



(12) 



Now lot us transform to magnetic normal coordinates 

fSAq,l)\ 
5.(9,2) 



S,{q, 1 
V 56(9, 2)7 



f^+,i;(g)S'+,.Jg) + cr_..(q)S_.^(q) 



- a+ ,(q)5+.b(9) + a_.b(g)5_.b(9). (13) 



Here, cr+,.(g), cr_ ,,(q), CT+^b(q) and cr__f,(q') are the mag- 
netic order parameters for a magnetic structure with 
wave vector q. The diagonal form of the magnetic free 
energy ([6|) is therefore 

fmaa = lY.l-^+Al^T)\a+,,{q)\^ + C_Jq,T)\a_Jq)\ 

+ c+,6(g,T) \<j^Aq)\^ + C,b(9,r) |^_,,(g)n 

+ 0(a4). (14) 

At high enough temperatures, the eigenvalues © are all 
positive and therefore the stable phase is the paramag- 
netic one. As we lower the temperature, we reach a criti- 
cal temperature for which one of the eigenvalues vanishes. 
We denote the wave vector for which one of the eigen- 
values vanishes first as qj^- Since A^(q) > A_(q) and 



D > 0, the first eigenvalue which reaches zero is C+ x- 

At the temperature Tj^^g at which C+ x = there is a 
phase transition from the paramagnetic phase to the AF3 
phase, in which a_^. ^{lic) 7^ ^ but all other order param- 
eters remain zero. At the second transition AF3— )-AF2, 
the order parameter bilic) oi'dcrs as well. This is 
true provided that 



(15) 



The last condition ensures that bilic^ vanishes be- 
fore C_ x(9/C'-^) ^^'^ temperature is lowered. Hence- 
forth, we will omit the plus sign in the order parameters 
subscript. 

To describe the electric polarization, we need to add an 
electric free energy and a magnetoelectric coupling term 
to the magnetic free energy. Assuming a homogeneous 
polarization, the expression for the electric free energy to 
lowest order is 
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cell / J 



(16) 



where "l^gj; is the volume of the unit cell, P is the 
ferroelectric order parameter and Xe a ^be high- 
temperature electric susceptibility along the a direction. 
By symmetry considerations,— the allowed magnetoelec- 
tric coupling term of the lowest order in the incommen- 
surate phases is 



k(,(g/c)l sin(¥'x - ^b)Ph 



(17) 



where (^j. and Lp^^ are the phases oi a ^.{(Iiq) and crj(g/c), 
respectively, and r is a small real magnetoelectric cou- 
pling parameter. Below we examine the results of the 
model. 



III. PHASE BOUNDARIES AND ORDER 
PARAMETERS 

A. MnW04 without magnetic fields 

The wave vector q^^ that characterizes the AF3 and 
AF2 phases is determined by maximizing the eigenvalue 
A^(q) for a given set of coupling energies {Jj}. After 
carrying out the maximization procedure, we can find 
the first transition temperature by equating C^^ ^ to zero 
for q = qic: 



(0) 



A+(q/c) + ^ 



7V3 



(18) 



The index indicates that this is the transition tempera- 
ture in the absence of external magnetic fields. By trans- 
forming to normal magnetic coordinates, the free energy 
of the incommensurate phases up to the fourth order in 
the magnetic order parameters is 



5 



f={aT-D- X+{qjc)) W.AQic)f + 36 k.(9/c)l' + {^T - X^Qjc)) \<ybiqic)f + 36 ^^(q/c)!' 

3 p 2 
„=1 ^AB^a 



(19) 



r 



This expression is obtained by keeping the Fourier com- 
ponents q = ±qic in the total free energy / — f^^g + 
fei + fint- Minimizing with respect to the polarization 
components, we find the induced polarization 



Inserting Eqs. (120]) into Eq. (|19|), we get 



CT:,(97c)l KiQic)\^M'Px ~ "Pb)- (20) 



cell 



f^{aT-D- X^qjc)) WAqic)\' + 36 KiQict + K - Kiqic)) WbiQic)f + 36 Wb{qic)\' 
+ 26 WAQicf WbiQic)f [2 + cos(2¥., - 2^J - 27sin2(^^ - ^J] , (21) 



where 7 is a dimcnsionlcss parameter given by 



(22) 



In order to minimize the free energy (PT|) . the phase difference (p^ — (/?^ should be ±7r/2. In addition, we show below 
that 7 is of order 10~^. Hence the last factor in the square brackets of Eq. (PT|) will be neglected in the description 
of the magnetic phase transitions. The minimization of the free energy ([2T|) with respect to \(^xi^ic)\ ^^'^ Wbi^ic)\ 
yields 



k°(9/c)| 



k°(9/c)| 



66 



yb{qjc)\=0 , T^;;l<T<T'^l, 



n(0) 



'a 4T, 



(0) 



T 



(0) 
N2 



3T) 



246 



K('^.c)h\/^ 



T < T; 



(0) 



JV2' 



(23) 



and the corresponding free energies are 



f(0) 
J AF?, 



AO) 

JAF2 



T) 



12b 

a2[4(T^")-T) 



^N2 ^ ^ ^ ^N3^ 



(0) 
N3 



'-N2) 



■(T 



N3 



T, 



(0)N 



N2) 



32b 



T <T, 



(0) 



W2' 



(24) 



with the transition temperature T^^\ given by 

"2a 

By calculating the phase (f>{qj(j) of J{qj(j]Ti, 
can find the magnetic structure of the phases AF3 and 
AF2. Using the experimental incommensurate wave vec- 
tor q^(j = (—0.214,0.5,0.457), this phase is found to be 



(0) 



N2 



rpiO) _ — 



(25) 
T2) we 



4'{Qic) ^ 27ry for the two sets of exchange couplings. Us- 
ing this relation and ip^ — ipj^ = ±7r/2 in Eqs. pl| and 
dH), the spins of the two Mn'^+ ions in the AF3 and AF2 
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phases are 

S{R + T,) = V2\a°iqjc)\cos {qic'R + i^)^ 
TV2|ab°(q,c)|sin(g,c-« + V')b, (26) 
S{R + T^) = -^/2\al{qic)\cos{qic • -R + ?A + A(/))x 
±\/2|CT5?(q^(^)|sin(q^c • JR + t/' + A0)b. (27) 

Here tp is an arbitrary phase and A0 = • (t2 — Ti) + 
4>{.Qic) ^ = '^lic with g/c; c being the c component 
of qjc- Using the experimental value qjc ^ = 0.457,— 
this phase is A0 = 0.45 77r. This is exactly the magnetic 
structure observed in neutron scattering studiesJ^ We 
emphasize that while group theoretical analysis yields 
several magnetic structures consistent with the crystal 
symmetries, the magnetic structure described by Eqs. 
P6p is the actual structure observed in experiments. 
The two possible signs correspond to the phase differ- 
ence if^ — ipf^ = ±7r/2 and represent spirals with opposite 
chirality. 



S{R + T^) X SiR + T^) = 
±2|ff2(g/c)|k°(g/c)|sin (A</))z. 



(28) 



Here z is a unit vector perpendicular to the spiral plane. 
Various studies reveal that the spin chirality is strongly 
correlated with the electric polarization and can be con- 
trolled by poling the polarization with an external electric 
field.— 1^ This observation is in agreement with the form 
(|20p of the electric polarization, in which tp^ — ipf^ changes 
sign together with P. 

Taking into account the magnetoelectric coupling in 
the description of the magnetic phase transitions will in- 
troduce small corrections to the transition temperature 
T^2 to the order parameters in the AF2 phase. As 
mentioned above, these corrections are governed by the 
dimensionless parameter 7 [see Eq. (P^ j. Using these 
corrections to the first order in 7. wc find that the elec- 
tric susceptibility takes the form 



XE,b 

x'eA 
X%h 



t>t: 



(0) 



p(0) 



n(0) ^ 



^3 



f ^°2^ < T < T; 

T^°l <T<T, 



f;(o) 
-g{T)^ 



,(0) 
JV2 



(0) 
N3 

(0) 
N2 



T <T 



(0) 



where t'^I is the shifted transition temperature: 



T(0) _ r^(O) 

^N2 ~ -^N2 



The function g{T) is 
9iT) 



f2iT) 



(29) 



(30) 



(31) 



where 



with Vi 



(0) 



^7.(0) 
3 ^3 3 



(0) 



^2' 



(32) 



737lr°2 ) and = 16 ( T^g 



.(0) 



97^(0)^(0) 85 rp{Q) 
e N2-^ N3 ^ 6 \^N2 

The first order phase transition AF2^'AF1 can be 
treated in the following way. Since the API phase is 
characterized by the commensurate wave vectors qc 1.2 = 

(±i,i,i), we calculate the free energy ffpi for this 
phase and then look for a temperature below which 

fAFi < /af2- Since 9^2 = -Qci + (0,1.1) wc need 
to consider only the Fourier components q ~ ^qc — 
±(i, i, i) in Eq. After some algebra we find the 

free energy 

f^{aT-D-X^{qc))\a,{qc)f 
+ b\a,{qct [3 + cos{4^~ Any)]. (33) 

Here (p is the phase of a^{q^), determined to be Tr(y + 
j) in order to minimize the free energy. Therefore the 
equilibrium order parameter and the corresponding free 
energy are 



KiQc)\ 



\+{qc) + D-aT 
46 



,(0) _ {K{qc) + D-aTY 



J AFl ^ ' 



(34) 



(35) 



For the commensurate wave vector q^ — (;|,^,^) we 
find the phase 4>{qc) = 27ry — tt of J{qc] Tj) for both 
sets of exchange couplings. Using this relation and Lp = 
7r(y -I- i) in Eqs. (|lT| and ©, the spins of the two Mn2+ 
ions in the AFl phase are 

S{R + Ti) = V2\al{qc)\ cos {qc-R+ J)x, 



SiR 



-V2\cTl{qc)\ cos {qc R 



(36) 



Equations p6|) describe a magnetic structure of the type 
ttii along both the a and c axes, in agreement with 
the structure observed in experiments^^ We note again 
that this is the observed structure out of the two possible 
structures suggested by group theory. 



The solution of the inequality f^pi < /Xf2 '-'^ ^^'^ 
form T < T^°l provided that 



f(0) 



e > max < 2 



I (1 - 7?) , ^ (1 - VW^) 



(37) 



where e 



D 



and rj 



transition temperature t'^\ is given by 



In this case, the 



T 



(0) 



-1 



4e - 3e^ 



4(2 (r]-l) + e) 



(38) 
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We study below the effects of magnetic field on the tran- 
sition temperatures. 



In the incommensurate phases AF3 and AF2, Eq. (|40)) 
is replaced by 



B. The effect of an external magnetic field 

The formalism presented above can be generalized 
to take into account the effect of a uniform exter- 
nal magnetic field h. This can be accomplished by 
adding to the free energy the Zeeman term = 
5Mb J2rJ2t=t,,t^ + t) • /i, or, equivalently^ 



diMl^(lic)\'+d2Ml^^ic)\ 



where d^^ = = 126 and = d^^^ 
sponding form in the AFl phase is 



, (45) 



46. The corre- 



fz^^=9^^B E S{S),r)-h. (39) 



TV 



Minimizing the free energy with respect to 5*0,(0, r) at the 
paramagnetic phase, we find the response to the external 
magnetic field 



(a = x,6), (40) 



with the magnetic susceptibility following a Curie- Weiss 
law 



xjr) 



{g^^Bf 



aT-D^-2Y.\l,J^ 



(41) 



Comparing Eq. (|4ip with the general Curie- Weiss lai 



.33 



3^5 



(42) 



5jO,r) 



1 I -.H(Qc)r 



(46) 



with = 126 and = 46. The ferromagnetic Fourier 
component at q = couples to the incommensurate 
and commensurate wave vectors through the fourth or- 
der term in Eq. 1^. This coupling modifies the coeffi- 
cients of the free energy expansion and, consequently, the 
transition temperatures. In the presence of an external 
magnetic field, the first two transition temperatures are 
(to second order in the magnetic field) 



(0) 
JV3 



1 - 12. 



1-4- 



h = h^it, 



h = /ifjb, 



(47) 



we identify the parameter a introduced in the expansion 
of the entropy [see Eq. ©] as 



3kr 



J(J + 1) 



(43) 



For Mn^+ ions with J = S = 5/2 this parameter is ayi^ = 
0.343fc^. The Curie- Weiss temperature is related to the 
exchange couplings and the anisotropy energy by: 



J(J+ 1) 
3k o 
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(44) 



(0) 
N2 

(0) 
N2 



with K = 



1 

n(0) 



1 - 16k 



h = /ijb. 



(48) 



For an external magnetic 



field along the easy axis direction, the inequality which 
determines the stability range of the AFl phase is 



T <T, 



(0) 
Nl 



86 



T-3 



Kile) 



D 



18 T 



\+{qc) + D 



8(t^)-t)' 



2T. 



(0) 



Af3 



9t^B 



(49) 



8 



while for a magnetic field along the b direction it is 



2(A+(9c)-A+(g,c)) +^ 



X+{qc) + D _ 2 / (0) (0) 



2 T 



4- 9t(°) 
^ jV3 ^''Af2 



(50) 



Equations (|T7|) - (|50|) describe the T ~ H phase diagrams 
up to second order in h. 



C. The effect of doping 

We can gain insight on the effect of small concentra- 
tions of magnetic Fc^+ or non-magnetic Zn^+ and Mg^+ 
ions at the Mn sites in the following way. Assuming 
that the orbital angular momentum is quenched, we set 
J = 5 = 2 in Eq. and identify the parameter a [see 



Eq. ([2])] for the ¥e^+ ion as 
value, we get 



0.5kg. Using this 



aix) =aM„2;-t-aFe(l-a;), 



(51) 



where x is the Fe concentration. Since the exchange cou- 
plings of Fe-Fe and Fe-Mn pairs as well as the anisotropy 
energy for the Fe ion are not known, we assume a lin- 
ear dependence of the quantities A_|_(q'j(^), D and ?/ = 



^+^^'^' for small values of x: 



A+(g/c(^)'^) A+ ^^(0) + Cia;, 
D{x) = D{{))+C2X, 
■q{x) = r/(0) + Cgx. 



(52) 



We use the relations (|52p in order to modify the expres- 
sions (dH]), (PSj) and for the transition temperatures. 
Then, by expanding these expressions to first order in x 
and fitting to the slopes measured in experiments;^ we 
are able to extract the values of c-^ , C2 and C3 . We neglect 
any changes in the parameter b. 

For the case of the non-magnetic Zn^^ ion we set a-^^ = 
D^^ = as weU as jf"-^" = jf"-^" = 0, and find the 
x-dependence of the different parameters 



>^+iq,x) = A+(q)(l - xf, 
a{x) = a(0)(l-a;), 
D{x) = D{0){l-x), 
r]{x) = 77(0). 



(53) 



Using these relations the first two transition tempera- 
tures are given by 



(0) 
^3 



A+(g/c) , 

a 

A+(97c) 



(54) 



These results explain the linear decrease of Tjyg and of 
TjYj as a function of x observed in experiments<^i^ The 
treatment of the AF2— >-AFl transition is much more sub- 
tle and will be discussed below. We note that all the 
results above do not depend on the type of the non- 
magnetic ion. This is in agreement with the observed 
similarities of the transition temperatures in Zn^"*" and 
Mg2+ doping^.!^ 



IV. COMPARISON WITH EXPERIMENTS 

In this section we compare the results of the preceding 
section with different experimental observations and ex- 
amine the consistency of the phase diagrams with the ex- 
perimental sets of exchange couplings of Ehrenberg et al. 
and Ye et al.. The results of the preceding section can be 
used to fit the parameters of the model within the Landau 
theory. We use Eqs. (HH) and ([25]) with Oj^j^ = 0.343/0^ 
[see Eq. (j43p ] and the experimental transition tempera- 



tures Tj^g and order to extract the values of the 

parameters \j^{qi(j) and D for MnW04. 



.(0) 



Using the ex- 

13.5ii: and T^°l « 12.3- 12.7is:, 
these parameters are found to be A^(7jp) ~ 4.36— 4.45i4r 



perimental values T, 



(0) 



NZ 



and D 0.27 - Q.lSk^K. The ratio ?/ 
chosen to be 77 



is then 



+ ^(lic) 

0.97 - 0.98 in order to fit Eq. ^ to 



the experimental transition temperature Tj 



(0) 



7 - m. 

These values are consistent with the condition (|57|) . 

Next we use Eqs. (gS]) and (|47l)-(l50l), with the exper- 
imental Curie- Weiss temperature 6^ ~ ~ —75K^^^ 
and calculate the T — H phase diagrams by fitting the pa- 
rameter b. In order to get the best fit to the experimental 
phase diagram of Arkcnbout et al.^ the parameter b was 
chosen to be 0.095fc^/-r. Figure [1] shows the results. The 
calculated and the experimental phase diagrams are in 
good agreement. Discrepancies at low temperatures or 
at high fields are expected due to the finite expansion of 
the free energy, which is terminated at fourth order. 

The development of the magnetic order parame- 
ters with decreasing temperature has been studied by 
polarized-neutron diffractions.— Generally, the magnetic 
moment at site r belonging to the unit cell at the lattice 
point R can be written as 



m{R + t) = cos {q R + (j)^) x 
-I- mf, sin {q- R + 0^) b. 



(55) 
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FIG. 1: (Color online) Magnetoelectric phase diagrams of 
MnW04 with magnetic fields parallel to the (a) easy and (b) 
b axes. The solid lines are the calculated transition tempera- 
tures and the dots describe the experimental points of Arken- 
bout et al.^ The calculated phase diagrams were obtained 
by setting h = 0.095fcsiC. 



The cross-sections for polarized-neutron scattering, 
where the neutrons are polarized parallel and anti- 
parallel to the scattering vector, are given bj«^ 

\2 



(56) 



with Iq being a constant. Using Eqs. (^5)) and ([M]). 
we see that these cross-sections are proportional to 

iH{(lic)\ ± H{Qic)\f + kS(gc)|'- Then, from the 



second of Eqs. (|23l) . the magnetic order parameters in 
the AF2 phase can be written as 



k°(9/c)| 



7^(0) _ r(O) 



^N2 ^ 



(0) 



N2 



T) 



(57) 



Toledano et al^ assumed that [^^(qj^)! is fixed below 
T^2- According to the first of Eqs. ([57| . such an assump- 
tion is valid only for r^°] - T < | (t^^I - Tj°},) . At 



lower temperatures this assumption is inconsistent with 
the evolution of the observed integrated intensities re- 
ported in Ref. [13, which show that both |cr°(q'j(^)| and 
|'^b('3'/c)| continue to grow below Tjyj, with the ellip- 



ticity p 



approaching 1 (so that the 



spiral is almost circular) as the temperature decreases. 
Therefore, we preferred to use the explicit dependence 
of |cr2(<7jp)| on the temperature. Using Eqs. (|57p . the 
ellipticity below T^, can be written as 



1 



P = 



(58) 



where i 



4(T«, 



Since the difference T^^ - 



'-^N2 



3(^2-^) 

0.8K is very small in the case of MnW04, the ellip- 
ticity rapidly approaches 1 with decreasing tempera- 
ture in the spiral phase AF2. The small difference 
r^3 — fo'^ MnW04 is a consequence of the small 
single-ion anisotropy of Mn^+ ions. This should be com- 
pared with the case of TbMn03, for which T, 



and T, 



N3 



A2K 



N2 



27K. In this multiferroic, the ellipticity 



grows much more slowly with decreasing temperature 



36 



due to the large difference 



N3 



^N2 



15K, which 



is in turn a result of the larger single-ion anisotropy 
of Mn'^+ ions. In Fig. [5] we sketch the quantities 

(K(g/C7)| ± kfc('7/c)|)'+k2(9c)|' from Eqs. ^ and 
(1341) together with the experimental data points of Ref. 




8 9 10 11 12 13 14 15 16 

Temperature (K) 



FIG. 2: (Color online) The temperature dependence of 

i\^UQic)\ ± k"(qic)|)' + kS(qc)r- The red (thick) line 
corresponds to the -|- sign and the blue (thin) one to the mi- 
nus sign. The integrated intensities of t he p olarized-neutron 
diffraction (scaled by 0.002) from Ref. [Sy are the red dots 
and the blue empty squares. 



The development of the calculated order parameters 
is in a qualitative agreement with the temperature de- 
pendence of the integrated intensities. However, for 



(0) 



JV2 



T > I (tI°1 - rj^l ) in the AF2 phase, the quan- 



tity (|o'°(<7jp)| -f- |cr°(g7c)|) is linear in T, in contradic- 
tion with the temperature dependence of the integrated 
intensity, as can be seen in Fig. [21 A possible explanation 
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for this apparent discrepancy is related to fluctuations 
near the transitions, that are not taken into account by 
the mean-field Landau theory.— As pointed out in Ref. 
Issl . the transition P— >AF3 belongs to the universality 
class of the XY model, while the transition AF3— >AF2 
belongs to the Ising universality class. Hence we present 
in Fig. [3] the same quantities as in Fig. [21 but replacing 
the square roots of Eqs. ([57|) and (|34p by the critical ex- 
ponent P — 1/3, roughly appropriate for these two mod- 
els. As seen from the figure, these revised expressions arc 
in good agreement with the observed integrated intensi- 
ties. This behavior illustrates the possible importance of 
fluctuations in MnW04. Further consequences of fluctu- 
ations near the transitions will be discussed below in the 
context of the Ginzburg criterion. 
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FIG. 3: (Color online) The temperature dependence of 

(kS(<7/c)| ± k°(9/c)|)^ + kS(<7c)r "it'^ critical expo- 
nent /3 ~ 1/3. The red (thick) line corresponds to the + 
sign and the blue (thin) one to the minus sign. The inte- 
grated intensities of the polarized-neutron diffraction (scaled 
by 0.0014) from Ref. [s^ are the red dots and the blue empty 
squares. 



FIG. 4; (Color online) (a) The ferroelectric polarization and 
(b) the dielectric constant The solid lines are the calcu- 
lated quantities and the dots are the data points of Taniguchi 
et a/..— The calculated polarization was obtained by setting 

4^ =21^iC/m2. 



The magnetoelectric coupling r is determined by fit- 
ting Eq. ([20)1 to the experimental data of the induced 
ferroelectric polarization The ferroelectric polarization 
is plotted in Fig. 4(a) The best fit to the experimental 

data is obtained for the value ^^■'•^'^^ = 21^C/rn?. In 

addition, the electric susceptibility for T > T^l (in the 
paraelectric and paramagnetic phase), is experimentally 
found to be x'e b ~ ll-Seg'^ The dimensionless param- 
eter 7 [see Eq. ([22])] is then 7 = 5.9 • 10~^. This value 
supports the assumption that the magnetic transitions 
are almost unaffected by the magnetoelectric coupling. 

x° 

The dielectric constant e,, = 1 -I — — is shown in Fig. 

^0 

[4(b)| This result is in good a gree ment with the experi- 
mental measurements of Ref. llSl The narrow width of 
the divergence region is a consequence of the small differ- 
ence between Tj^j xj^l . Once again, the discrepancy 
between the linear behavior of the calculated polarization 
and the observed one may be reconciled by assuming a 
critical exponent /3 ~ ^ for the magnetic order parame- 



ters. The behavior of the calculated polarization in this 
case is given in Fig. [5] 

To examine the effect of Fc doping, we use the rela- 
tions ([?T|) and ([5^ in the expressions for the transition 
temperatures and fit the slope to the experimental value 
according to the x — T phase diagram of Chaudhury et 
al.^ This procedure yields the values « —3.26kgK , 
C2 ~ IS.OSkgK and C3 ~ —1.3. The anisotropy energy 
increases with increasing Fe concentration, as expected, 
since as opposed to the Mn^+ ion, the Fe^^" ion possesses 
a non-vanishing angular momentum)^ 

Calculating the different parameters for a small Fe con- 
centration X and repeating the calculations of the T ^ H 
phase diagram, we can check the consistency of the above 
results. The resulting phase diagram for x = 0.035 is 
shown in Fig. [51 Except for high fields or low tem- 
peratures, the result is in fine agreement with the mea- 
surement of Ye et al.^ The reentrant ferroelectric phase 
observed at low temperature o^°i"^° may be explained by 
higher order terms in the free energy expansion. 

The effect of non-magnetic ions on the transition tcm- 
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FIG. 5: (Color online) The ferroelectric polarization calcu- 
lated with the critical exponent /? ~ i. The solid lines are 
the calculated quantities and the dots are the experimental 
points of Taniguchi et aZ..— The calculated polarization was 

obtained by setting y-'' = 27.5fiC/m^. 
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FIG. 6: (Color online) Magnetoelectric phase diagram of 
Mno.965Feo.o35W04 with the magnetic field parallel to the 
easy axis. The solid lines are the calculated transition tem- 
peratures and the dots are the experimental points of Ye et 



peratures Tpj^{x) and 7\r2(^) given by Eq. ([54)) . These 
results are drawn in Fig. [7] together with the experimen- 
tal data of Chaudhury et al^ of Mnj_^Zn^W04. Simi- 
lar results have been observed in Mn;^_^Mg^W04<^ We 
stress that unlike the case of Fe doping, the results for the 
transition temperatures Tj^^(x) and Tj^2{^) i^i the case 
of non-magnetic ions doping do not require additional 
phenomenological parameters. 

As opposed to T^g (x) and {^) ? the calculated tran- 
sition temperature T^j^(x) does not coincide with the 
experimentally measured onci^ The discrepancy may 
be explained by allowing small changes in the exchange 
couplings J^^n-Mn spin-latticc coupling (or ex- 

change striction). In other words, if we assume that 
jM„-Mn(^) = Jf^'^-^'^il + i^x) with < 1, then T^^{x) 
changes dramatically while T^^{x) and Tj^^i^) 
most not influenced. The reason for this behavior is that 
the transition temperature Tj^-^ [see Eq. psp ] is much 



more sensitive to small changes in the exchange couplings 
than the transition temperatures Tjyg and Tpf^ [see Eqs. 
([TS]) and (pSj)]. A significant spin-lattice coupling in the 
multiferroic MnW04 has been demonstrated^ by the ap- 
pearance of an incommensurate lattice modulation in the 
AF3 and AF2 phases, with a lattice propagation vec- 
tor equal to twice the magnetic propagation vector. In 
addition, thermal expansion measurements reveal con- 
siderable discontinuities in the lattice parameters at the 
AF2— >-AFl first order phase transitions^ Another indica- 
tion for a dependence of the Mn-Mn exchange couplings 
on the non-magnetic dopant concentration is provided 
by the small change of the incommensurate propaga- 
tion vector from Qj^j (-0.214,0.5,0.457) in MnW04 
to ^ (-0.209,0.5,0.453) in Mug ggZuQ i5W04^ 
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0.1 0.2 0.3 0.4 0.5 

Zn Concentration (x) 

FIG. 7: (Color online) Transition temperatures Tff^{x) and 
Tn2{x) of Mni_^Zn^W04. The solid lines are the calculated 
transition temperatures and the dots are the data points of 
Chaudhury et al.^ 

The next step is to compare the above fitted param- 
eters with the parameters calculated directly from the 
experimental sets of exchange couplings of Ehrenberg et 
al. and Ye et ai. The calculated exchange couplings of 
Ref. [13 yield much higher transition temperatures than 
the observed ones and thus will not be discussed here. 
Indeed, The problem of overestimation of exchange in- 
teractions by DFT calculations has been indicated by 
the authorsj^i 

The first step is to maximize A^(q) [see Eq. (|lll) ] 
in order to find the incommensurate wave vector 
and the corresponding eigenvalue A_|_(qjp). The max- 
imization process yields Qjq ~ (—0.28,0.5,0.44) and 
A^(<7jp) = 3.82kgK for the — Jg values of Ehrenberg 
et al.^ while for the Ji — Jn values of Ye et al^ we find 
= (-0.3,0.5,0.49) and X+{qjc) = SMk^K. These 
results are in qualitative agreement with the incommen- 
surate wave vector = (—0.214,^,0.457) observed in 
experiments. However, the differences are not negligi- 
ble, suggesting possible errors in the experimental sets 
of exchange couplings. In addition, the transition tem- 
peratures Tjl^g and Tj:^l calculated from Eqs. (fT8| and 
(HH) with = 0.343fc^ [see Eq. (gS])] are found to be 

rj^l = 12. 79/^, T^^l = lO.SK for the set of Ehrenberg et 
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al. and T, 



(0) 



13.67K, T; 



(0) 



10 K for the set of Ye et 



N3 ~ it , J-JY2 

ai. These values shghtly differ from the observed transi- 
tion temperatures, especially the second one. The ratio 

V = t't^ is found to be 77 = 0.974 and 77 ^ 0.97 for the 
ylic I 

sets of Ehrcnbcrg et al. and Yc et ai, respectively. Table 
im summarizes the values of X^Iqjq), D and 77 calculated 
from the experimental sets of magnetic parameters and 
those fitted to the experimental transition temperatures. 



TABLE 11: Comparison between the model parameters calcu- 
lated from the experimental sets of Ref. [2y and Ref. [2^ and 
those fitted to the experimental transition temperatures. 



Parameter 


Ref. 26 


Ref. 28 


This work 




3.82 


3.85 


4.36-4.45 


D{ksK) 


0.568 


0.83 


0.27-0.18 


V 


0.974 


0.97 


0.97-0.98 



The calculation of the Curie- Weiss temperature re- 
veals a much more serious discrepancy. According to 
Eq. p4|) . the Curie- Weiss temperature is 0^ — — 7.6A', 
= for the set of Ehrenberg et al. and 



-23.2i^, 



-25.65K for the set of Ye et 



These values do not fit the experimental Curie- Weiss 
temperature 9 « —75KM^ We suspect that the origin 
of most of the discrepancies are errors in the set of mag- 
netic couplings. The results suggested by our model may 
be used as additional constraints in the determination 
of those couplings. As mentioned before, an additional 
possible cause for the above discrepancies is related to 
fluctuations near the transitions, as will be discussed in 
the next section. 



A^(q) to second order around qj(j: 
X+iq) ~ X+iqjc) + M,^ {q^ - qjc^, 



QIC, 



withAf„. = 



(60) 

Denoting by /.jj and /ig 



the three eigenvalues of the positive matrix — A/^-, the 
three correlation lengths are 



Mi 



(61) 



Substituting q = qj(. and N = in Eq. §9^, the 

condition {\Sa^{qj(j)\^) ^ |f°(9/c)r validity of 

the mean-field theory reads^i 



(0) 



B-^N3 



cell 



4a M-T) ^1^3 



6b 



(62) 



Inserting Eq. ((6T|) into Eq. ((62)) at the Ginzburg tem- 
perature Tq, wc find 



rp rp{0) 

^N3 



(63) 



Equation (j63p estimates the temperature range below 
Tjl^g, in which fluctuations are not negligible. 

Let us now estimate the Ginzburg range according to 
the second method. On the one hand, according to Lan- 
dau theory, the heat capacity c = — grows discon- 
tinuously at the transition P^AF3: 



V. THE GINZBURG CRITERION 

The results of the preceding sections have been ob- 
tained within the mean-field approximation. Here we 
estimate the Ginzburg range, in which fluctuations be- 
come important, near the flrst transition P— >AF3, by 
two methods. First we compare the mean square fluctu- 
ation of the order parameter (^xiQic) '^ith the mean-field 
value, and then we compare the discontinuity in the heat 
capacity derived from the Landau theory with the diver- 
gent heat capacity, originating from the fluctuations at 
quadratic order 1^ 

Let us denote by Sa^{q) = cr^{q) — {a^{q)) the fluc- 
tuation of the order parameter in the AF3 phase. The 
correlation function of these deviations is 



{6a,{q)S<yAQ')) 



^s'^^q' -q 



AN [D + \^{q)-aT) 



(59) 



where N is the number of unit cells in the correlation vol- 
ume. We can find the correlation lengths by expanding 



^N3 

6b ' 



(64) 



On the other hand, assuming fluctuations at quadratic 
order, the singular part of the heat capacity is given by 



2(27r)^ 



BZ {aT-D- X^{q)y 



(65) 



where the integral is over the first Brillouin zone. In 
the neighborhood of , the main contribution to the 
integral comes from the neighborhood of the incommen- 
surate wave vector qj(-, in reciprocal space. Thus we can 
use the expansion (|60p . Replacing the first Brillouin zone 

by a sphere, and taking T w T^g, we can estimate the 
integral in Eq. (|65p : 



-0.5 



(66) 
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Comparing Eqs. (|64l) and (|66l) at the Ginzburg temper- 



ature Tq, we find^i 



T 



f 6 y fc|^'Ke» 

V167r j 



(67) 



Calculating the eigenvalues ^j^, /ij ^'^'^ Ms froni the ex- 
perimental sets of exchange couplings, the Ginzburg tem- 



perature is estimated to be 



(0) 



JV3 



9AIK and 



Tn-T, 



(0) 



N3 



6.24:K for the sets of Ehrenberg et al. and 
Ye et al, respectively, by the first method [see Eq. ([M])] 



while it is 



(0) 



N?, 



0.04A' 



i 0.06A' and 

by the second method [see Eq. (|67|)]. These values sug- 
gest that fluctuations of the order parameters can also 
contribute to the discrepancies between the experimen- 
tal data and the mean-field Landau theory results. 



VI. SUMMARY AND CONCLUSIONS 

We have studied the phase diagram of Mi\i_.J>A^WO^ 
(M=Fe, Zn, Mg) by a semi-phenomenological Landau 
theory. The energy has been modelled by a Heisenberg 
Hamiltonian with a single-ion anisotropy, while the en- 
tropy has been expanded in powers of the classical spins. 
This approach is different from the previous theoretical 
studies^i^ which are purely phenomenological, since it 
enables to compare different sets of exchange couplings. 
Although a purely phenomenological approach may cap- 
ture all the symmetry aspects of the problem and may 
provide a full mapping of the stable states allowed by 
the order parameter symmetries,— it does not indicate 
a clear connection between the free energy coefficients 
and the microscopic interactions. The advantage of our 
approach is the simple relation of the free energy coeffi- 
cients with experimentally derived quantities such as the 
supcrcxchange couplings and the anisotropy coefficients. 
For instance, this simple relation allows us to consider 
the effect of different dopants on the phase diagram, not 
discussed in Ref. [H. We emphasize that our approach 
does not contradict any symmetry requirement. 

We used the superexchange interaction couplings from 



the inelastic neutron scattering studies of Ehrenberg et 
al^ and Ye et al.^ The results show that both sets 
yield transition temperatures T^g and Tj^l that slightly 
deviate from the experimental temperatures, and signif- 
icantly underestimate the Curie- Weiss temperature \9\. 
In addition, the calculated incommensurate wave vector 
qj(j has non-negligible deviations from the experimen- 
tally observed one. The results presented here can serve 
as additional constraints on a future determination of 
the magnetic Hamiltonian parameters. Another possible 
cause for the discrepancies relates to fluctuations near the 
transitions. We have demonstrated the possible impor- 
tant contribution of fluctuations in MnW04. This issue 
should be further examined in future experiments. 

Beyond that, the model clarifies the effect of different 
dopants on the phase diagram. The sensitivity of the 
expression (|38p for the transition temperature Tj^^ix) 

to small changes of the ratio ri = reflects the 

frustrated nature of the multiferroic MnW04. The ori- 
gin of the complex phase diagram lies in the competi- 
tion between different superexchange interactions. Small 
changes in the local environment of the Mn^"*" ions due to 
a chemical doping cause a signiflcant change in the phase 
diagram. The sensitivity for the local environment man- 
ifests itself by the contrasting behavior of doping with 
different ions. 

Looking to the future, two points should be further 
examined. Firstly, a new analysis of the inelastic scatter- 
ing experiments, together with the additional constraints 
provided in this work, should improve the exchange cou- 
plings for the multiferroic MnW04. Secondly, the mea- 
surement of the critical exponents near the transitions 
would shed light on the effect of fluctuations. This may 
contribute to the general understanding of critical phe- 
nomena in multiferroics. 
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